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INFLUENCE OF OIL-SQUEEZE-FILM DAMPING ON STEADY-STATE RESPONSE 
OF FLEXIBLE ROTOR OPERATING TO SUPERCRITICAL SPEEDS 
by Robert E. Cunningham 
Lewis Research Center 

SUMMARY 

The response of a long, slender, ball-bearing-supported rotor was measured for 
unbalances ranging from 0. 62 to 15. 1 gram -centimeters. Three unbalance distributions 
were also evaluated. The test apparatus contained squeeze-film dampers at the bear- 
ings, and the rotor could be operated with or without the dampers being activated. Re- 
sults of undamped and damped responses were compared with an unbalance response 
computer program by using a mathematical model which closely simulated the actual 
hardware. 

Squeeze-film damping coefficients were calculated from measured data and com- 
pared with theoretical short -journal -bearing approximations over a frequency range from 
5000 to 31 000 cycles per minute. The dampers were fed at the center from a circum- 
ferential oil distribution groove and had a length -diameter ratio of 0, 094. 

Experimental dampii^ coefficients were in good agreement with values calculated 
from short -journal -bearir^ theory. The values generally remained constant over the 
frequency range investigated. 

Squeeze-film dampers located at the bearings were very effective in reducing am- 
plitudes of motion at the rotor ends. As much as 88 percent reduction in amplitude was 
achieved at resonance for a total imbalance of 0. 93 gram -centimeter distributed in a line 
and in the same phase over three disks. 

When oil -squeeze-film dampers were activated for an in-line, alternating -phase un- 
balance distribution, three distinct resonant amplitudes were reduced to only one. The 
mode shape for this one corresponded to the third bending critical. 

The oil-squeeze-film dampers were also effective in reducing the amplitudes of mo- 
tion for nonsjuichronous as well as synchronous whirl. 

An unbalance response computer program was quite accurate in predicting the re- 
sponse of the experimental rotor having an in-line, in-phase unbalance distribution over 
the three disks. 

Operating the rotor at large eccentricities induced a damper instability character- 
ized by a sudden increase in the eccentricity ratio. This phenomenon has been observed 
by other investigators and is referred to in the literature as bistable operation. 


INTRODUCTION 


Continuing efforts to improve the performance of aircraft turbojet engines are re- 
sulting in higher and higher operating speeds. Such applications as small helicopters, 
remotely piloted vehicles (RPV’s) and small private and commercial aircraft are calling 
for increased specific -thrust ratios by reducing the overall weight of powerplants 
(refs. 1 and 2). Aircraft engine manufacturers are also looking at ways of reducing the 
power losses in an engine. One obvious loss is the leakage of air between compressor 
stages. Such losses could be reduced by decreasing the seal clearances, particularly 
those at compressor blade tips. 

Higher rotor operatii^ speeds, lighter powerplant structures, and smaller blade tip 
clearances will call for greater emphasis on the dynamics of rotors and their bearii^s 
and on the maimer in which the bearings are mounted in the frame. The combination of 
lighter rotors and higher operating speeds can result in operation above one or more lat- 
eral bending critical speeds. Vibration amplitudes at these critical frequencies can be 
quite high. The reduction of these vibration amplitudes to an acceptable level will dictate 
the use of more sophisticated balancing methods. One method currently receiving a 
great amount of attention is the high-speed, multiplane balancing method described in 
reference 3. 

Sometimes, however, it is not possible to achieve low levels of vibration by bal- 
ancii^ alone, particularly with lightly damped ball bearings and the built-up rotor as- 
semblies which compose the modern turbojet engine. These rotors, thoT^h balanced 
well initially, degrade with continued service time. Compressor and fan blades can be 
damaged by debris ingested at the intake. Loosening of tie bolts in built-up rotor assem- 
blies can chaise the stiffness of the rotor. Normal wear of rotating components also 
contributes to a gradual degradation in the balance of the rotors. Thus, techniques other 
than balancing need to be employed to control both steady -state and transient vibration 
amplitudes. 

Coulomb (dry friction) and viscous -oil shear mechanical dampers have been used for 
many years to control vibration amplitudes. Normally they have been used in controllii^ 
torsional vibrations in all types of rotating machinery (ref. 4). The inherent damping 
properties of hydrodynamic oil-lxibricated bearings were investigated by Ha^ and Sankey 
in 1956 (ref. 5). 

It is a well known fact that ball and roller bearings contain very little inherent damp- 
ing, and such bearings are used almost exclusively in aircraft turbine engines. The 
availability of lubricating oils and the requirement of large amounts of energy dissipation 
in a limited space were natural inducements to the ultimate development and use of oil- 
squeeze-film dampers (ref. 6). 

To properly design such a damper, the machinery design engineer must know some- 
thing about the dynamic characteristics of the rotating system. Normally he has a 
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niunber of computer programs available to him. Such pr(%rams can provide him with 
the flexible critical speeds and mode shapes of a particular rotor model in graphical out- 
put (e. g. , ref, 7), The computer program of reference 7 will also draw, to scale, a 
cross section of the rotor under study and will plot critical speeds as a function of sup- 
port stiffness. Also available are unbalance response computer programs that can pre- 
dict the response of almost any type of rotor on flexible, damped supports. Damping 
coefficients can be an input to the program either as constants or as polynomial functions 
of rotor speed (ref, 8), The output of such programs, however, is only as accurate as 
the values of input (e, g, , support stiffness and damping characteristics), 

A number of anal 3 d:ical papers have been written which develop the equations for cal- 
culating squeeze-film damping and stifhiess coefficient; two such papers using the short- 
bearing approximation are references 9 and 10, One of the objectives of this investiga- 
tion was to determine experimentally damping coefficients for an oil -squeeze -film 
damper and to compare them with theory, 

A second objective of this investigation was to examine the experimental response of 
a three -disk flexible rotor to varying amounts of unbalance with and without squeeze -film 
damping at the bearing supports. The resvilts would be compared with an imbalance - 
response computer program similar to the one described in reference 8, This program, 
however, was modified to consider only circular orbits, and an iterative solution was 
added that includes support damping as a function of journal eccentricity ratio. 

Unbalances tested varied from a low of 0. 62 gram -centimeter to a hi^ of 15. 1 
gram -centimeters. Three different types of unbalance distributions were investigated: 

(1) T37pe I. Identical set screws were placed in the same radial plane in each of the 
three rotor disks. The set screw in the center disk, however, was out of phase by 180°. 

(2) Type n. Identical set screws were placed in the same radial plane in each of the 
two end disks only. The two set screws were 180° out of phase. 

(3) T 3 q)e ICE. Identical set screws were placed in the same radial plane in all three 
disks. Results were obtained for a rotor speed range from 5000 to 31 000 rpm. 


APPARATUS 
Mechanical Features 

Experimental data were obtained with the apparatus shown in figures 1 to 5. The ap- 
paratus consists of a long slender shaft supported in ball bearings which are 48. 3 centi- 
meters apart. A single, integral di^ is located at midspan, while removable disks of 
identical size, one at each end (disks 1 and 3 in fig. 2), are located outside the bearings. 
The rotor is 69. 09 centimeters long and has a uniform diameter of 2. 54 centimeters be- 
tween bearings. The disks are 9. 5 centimeters in diameter and contain threaded holes 
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Figure 1. - Test apparatus used in experiments on steady-state unbalance response. 


every 10° on a ball circle of 8. 26 centimeters. Set screws can be placed in the holes as 
correction weights for balancing the rotor or to deliberately imbalance the rotor in a par- 
ticular bending mode. 

It was anticipated that severe bending strains would be induced in the rotor, and 
therefore it is made of a high-strength alloy steel, AISI 4340, which has a minimum yield 
stress of 1585 meganewtons per square meter. 

The rotor is supported in ball bearings, one of which is shown in figure 4. The outer 
races of these bearings are located in cylindrical steel housings, which also serve as 
damper journals (figs. 3 and 4). Also shown in figures 3 and 4 are the aluminum center- 
ing springs. They provide flexibility in the bearing support, prevent rotation of the 
damper journal, and center the journal in the damper bearii^. The aluminum used for 
these springs is 6061-T6, which has a yield stress of 296 meganewtons per square meter. 

In order to determine the lateral stiffness, the springs were installed in the test rig 
with the dampers attached and then subjected to known loads in two planes. Deflections 
were measured and stiffness was calculated from the load, deflection curves. Slopes of 
these curves were linear, and the average value of the stiffness was 2. 36 meganewtons 
per meter. Values of stiffness in each of the two planes differed by only one-half of 
1 percent. 

The damper used in this investigation is similar in design to the one described in 
reference 11, with the followirg exceptions: No piston ring seals were used at the 
damper ends, although the grooves were retained (see fig. 3). The piston ring seals 
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Figure 2 . - Schematic of test apparatus used in experiments on steady-state unbaiance response. 


were eliminated because they introduced a tare damping and stiffness that could not be 
readily quantified. Eliminating the end seals, however, necessitated a reduction in 
damper radial clearance from 127 micrometers to a nominal clearance of 64 micro- 
meters in order to keep damper oil flows within the pumping capacity of the oil supply 
system. 

The radial, deep-groove ball bearings are 204 series bearings, having a 20- 
millimeter bore. They are designated high-speed turbine bearings of aircraft quality 
and have lightwe^ht phenolic retainers. The outer race fits slightly loosely in its hous- 
ing, and a preload is applied with wave washers. Air -oil mist lubrication for the ball 
bearings is supplied through three nozzles. The lubricant is a typical turbojet engine 
oil, type n ester, meeting the specification MIL -L -23 699A. The same oil is used in the 
squeeze-film dampers. 

The journals for the squeeze -film dampers are AISI 440-C steel hardened to 
Rockwell C 58 to 62. The jovirnal outside diameter is 7, 9 centimeters, and the width of 


5 



BALL BEARING 
HOUSING AND 
DAMPER JOURNAL-s 


20 mm BORE 
BALLBEARING 


DAMPING 
OIL IN 


r DAMPER 
CENTERING 
SPRING 



DAMPER i 
BEARING-' 


-SQUEEZE FILM 

'-CIRCUMFERENTIAL OIL 
SUPPLY GROOVE 


R. DAMPER JOURNAL RADIUS, 3.96 cm 
L, DAMPER JOURNAL LAND LENGTH, 0.74 cm 
Cr, DAMPER RADIAL CLfARANCE. 63. 5 pm 


Figure 3. -* Schematic of oil -squeeze-film damper used in experiments 
for steady-state response of three-disk rotor. 



Figure 4. - Components of flexible damped support. 


each of the two lands is 0. 74 centimeter, so that the length -diameter ratio is 0. 094. 

Oil is supplied to the two halves of the damper by a circumferential groove 0.64 centi- 
meter wide. The journal outer surface has a 0. 2-micrometer finish. The average 
assembled radial clearance in the dampers is 63. 5 micrometers. 

The outer damper bearing is held stationary (figs, 2 and 3). It is made of SAE 4340 
steel hardened to Rockwell C 45 to 50 and has a 0. 2 -micrometer finish. The flow of oil 
to each damper can be controlled independently with needle valves located in each of the 
supply lines. Temper ahme of the oil can be maintained within 2° C over the elapsed time 
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Figure 5. - End view of test apparatus showing ball bearing 2and oil- 
squeeze-film damper. 


of a typical test run. 

The test rotor is driven by a 10 -centimeter -diameter air turbine having a maximum 
operating speed of 60 000 rpm. The turbine is connected to the test rotor by a 0. 95- 
centimeter -diameter quill shaft and a spline coupling. 


Instrum entation 

Deflections of the test rotor from a neutral axis were sensed by noncontacting, eddy- 
current probes. These probes, aloi^ with signal conditioners and power siq)plies, pro- 
duced a voltage linearly proportional to distance. The approximate linear range for the 
probes used in this investigation was 900 micrometers. Before any of the eddy-current 
distance sensors were installed in the test apparatus, each was calibrated with a microm- 
eter fixture, and voltage oulput was read on a digital voltmeter. The calibrating fixture 
was designed by the probe manufacturer specifically for this purpose. A dynamic cal- 
ibration of several of the probes was made at one distance setting. The same calibrating 
fixture had a flat plate that could be rotated with a preset wobble underneath the probe be- 
ing calibrated. The output voltage was recorded and compared with the static calibration; 
agreement in all cases was good. 
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Sensors for monitoring shaft motions were located in mutually perpendicular axes at 
each of the three disks. Distance sensor pairs were also cemented into the damper bear- 
ing flush with its surface (figs. 2 and 6). These sensor pairs were displaced 45° from 
the damper force transducers. The force transducers, located in mutually perpendicular 
planes, can be seen in figure 5, 

The transducers used for measuring the forces transmitted through the damper oil 
squeeze film are quartz crystal, piezoelectric load washers. The load range is 0 to 
35. 6x10^ newtons with a stiffness of 175X10® newtons per meter and a resonant frequency 
of 50 000 hertz. Linearity over the load range was 1 percent, 

A fiber optic sensor was used to establish a reference signal for the measurement of 
phase angles, A piece of reflective tape was epoxied to the rotor surface with one edge 
coincident with the centerline of one of the balancing holes in the center dids (figs. 1 
and 2). The tape was 1, 3 centimeters wide and covered exactly 180° of the rotor surface. 
The square-wave voltage signal produced when the shaft was rotated was used as the ref- 
erence signal to which other sensor outputs were compared. 

A magnetic sensor near a six -toothed wheel moimted on the ttirbine quill shaft meas- 
ured rotational frequency (fig. 2), Ball-bearing temperatures were measured with iron- 

constantan thermocouples and recorded on a strip chart. Damper -oil -film temperatures 
were monitored on a digital readout. Oil flow rates to the squeeze-film dampers were 
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measured by turbine vane flowmeters located in the supply lines and were read otrt as 
frequency on a digital counter. 

Leads from the ntunerous sensors located on the test apparatus were connected to a 
patch board in the control room. From this board leads were paralleled to any number 
of readout and data acquisition systems. 

For example, the output signals from any one of the displacement or force trans- 
ducers could be directed into either of two channels of a trackii^ filter. The filtered out- 
put signal, a true sine wave at the rotor frequency, was paralleled to one input terminal 
of a phase meter and then to a peak-to-peak measuring voltmeter. In the phase meter it 
was compared with the square waveform generated by the fiber optic sensor. A phase 
angle was then detected and displayed digitally on a meter and also fed to the Y* axis of 
an X-Y-Y* plotter. The peak-to-peak voltage, representii^ either amplitude of motion 
or transmitted force, was displayed on a meter while simviltaneously being recorded on 
the Y axis of the plotter. Unfiltered signals could also be monitored with any of the 
readout instruments. Four X-Y curve tracing oscilloscopes were used to monitor si- 
multaneously any of the coordinate axis displacements or transmitted forces. 


DISCUSSION OF RESULTS 

The results of this investigation are presented in figures 7 to 16 and in table I, and 
are discussed in two sections. The first examines the response of a three-disk rotor to 
various magnitudes of unbalance. Three types of unbalance distribution were tested, and 
resiilts are shown in table I. The second section of the discussion compares ejq)eri- 
mentally determined damping coefficients with those predicted by theory for unbalances 
ranging from 0. 62 to 15. 1 gram -centimeters and rotor frequencies up to 31 000 cycles 
per minute. These results are plotted in figures 13 to 16. 


Unbalance Response 

Before any tests were conducted, it was necessary to determine the undamped crit- 
ical speeds and mode shapes. These were determined from the computer program of 
reference 7 and are plotted in figure 7 for a spring stiffness of 2, 36 meganewtons per 
meter. The deflections were only relative but indicated that oil -squeeze dampers placed 
at the supports should be effective when operatii^ through all fotir bending critical 
speeds. 

Before any deliberate unbalance distributions could be tried, the rotor had to be well 
balanced. First attempts to balance the rotor without any squeeze-film damping at the 
supports indicated a high d^ree of sensitivity of a ball -bearing -supported rotor to very 


9 


TABLE L - SUMMARY OF EXPERIMENTAL AND THEORETICAL DAMPING COEFFICIENTS 

AT VARYING TYPES AND MAGNITUDES OF UNBALANCE 

[Damper geometry shown in fig. 3-, type of oil, type II ester meeting specification 
MIL-L-23699-A; viscosity at 38° C, 0.033 (N) (sec)/m^; oil flow, 5. 15x10"“^ 

o 

m /min; supply pressure, 0.19 MPa.] 


Type of unbalance distribution 

Magnitude of 
total unbalance, 

u, 

g-cm 

Frequency range, 
cpm 

Average damping coefficient, 

®7)=0’ 

(N)(sec)/m 

Experimental 

Theoretical 

Type I - In-line, alternating- 

0.93 

11 000 - 31 000 

6740 

6492 

phase, three disks 

f 

1.67 

10 000 - 31 000 

6531 

6213 






Type II - In-line, alternating- 

0.62 

7 000 - 26 000 

6029 

6213 

phase, two end disks 





k 

1.11 

7 000 - 28 000 

6265 

5408 

1 

1.5 

7 000 - 28 000 

6843 

6212 


3.3 

6 000 - 29 000 

5416 

5880 






Type in - In-line, same-phase, 

0.93 

6 000 - 22 000 

5478 

57 58 

three disks 1 






1.67 

7 000 - 24 000 

5708 

5425 


2.25 

6 000 - 23 000 

5926 

5425 


4.97 

5 000 - 22 000 

5880 

5425 


10. 6 

6 000 - 30 000 

5289 

6108 


15.1 

5 000 - 25 000 

6375 

6843 
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Figure 7. - Undamped critical speeds and mode shapes. Centering spring stiff- 
ness, K, 23.63x10^ newtons per meter. 


small amounts of unbalance. Fairly good balance, however, was achieved by using a 
combination of the influence coefficient method (ref. 3) and single plane balancing. The 
results of rotor balancing, presented in figure 8, show that peak-to-peak amplitudes of 
motion at resonance did not exceed 25 micrometers at any of the three disks. A certain 
amount of machining runout and/or bend was present in the rotor, as shown by this plot. 
However, it was not possible to reduce this kind of runout further. 
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ROTOR FREQUENCY, cpm 

Figure 8. - Final balance of three-disk rotor. Oil dampers activated. 

Figures 9 to 13 compare the undamped rotor response with the response when oil- 
squeeze-film dampers were activated. The traces in figures 9 to 12 are for the filtered 
signal generated by either proximity distance sensors or piezoelectric crystal force 
transducers. The 10- and 20 -hertz -bandwidth filters track on the rotor frequency, and, 
therefore, the trace plotted is that of the rotor synchronous response only. Comparisons 
between undamped and damped response are for rather small unbalances of 0, 62 and 0. 93 
gram -centimeter because operating the rotor without any external damping at the bear - 
ii^s could have resulted in damaging the rotor and stationary seals if larger vihbalances 
had been iised. 

The plots of amplitude against frequency in figure 9 were produced by an in-line, 
alternating -phase distribtttion of weights in the three disks. The unbalance weights in the 
end di^s were placed 180° out of phase with the weight in the center disk, as depicted by 
the inset drawing. The total unbalance was 0. 93 gram -centimeter. The three set 
screws used as unbalance weights were carefully matched in weight and did not vary 
more than 4 milligrams. 

The amplitudes of motion plotted in figure 9(a) are those of the center disk. A max- 
imum peak-to-peak amplitude of approximately 230 micrometers was recorded for the 
rotor operatii^ with no externally applied damping. 
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ROTOR FREQUENCY, cpm 


(a) Center disk. 

Figure 9. - Comparison of rotor response for undamped and damped 
conditions. Type I unbalance distribution; unbalance, U=0.93 
gram-centimeter. 


The undamped rotor exhibited three distinct amplitude buildups at frequencies close 
to the ones predicted by the critical speeds program of reference 7, The maximum un- 
damped resonant amplitude occurred at 8100 cycles per minute. The type of unbalance 
distribution used in this trial wotild normally excite the first and third bending critical 
speeds, as can be seen from the mode shapes of figure 7. Therefore, it is reasonable to 
e:q)ect larger amplitudes to occur at resonant frequencies corresponding to those modes. 
It should be pointed out that, even though no external squeeze-film damping was present 
in this trial, it is obvious that some amount of damping is present in the oil -lubricated 
ball bearings. The maximmn amplitude predicted by an unbalance response program, 
for a ball bearing damping of 3500 newton -seconds per meter, is 381 micrometers at 
7500 cycles per minute. 

The lower portion of figure 9(a), for the test where oil-squeeze-film dampers were 
activated, shows one distinct buildup in rotor amplitude instead of three. The mode 
shape for this damped critical, occurring at 10 000 cycles per minute, was similar to the 
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0>) Disk 3. 


Figure 9. - Concluded. 


one at the third bending critical (fig. 7). The maximum peak-to-peak amplitude, when 
dampers were activated, was 216 micrometers, and this obviously was not a significant 
improvement. The tmbalance -response program, however, predicted a maximum am- 
plitude of only 85 micrometers at 9600 cycles per minute. The agreement between ex- 
periment and theory for this unbalance distribution was obviously not very good. Because 
of the large rotor exciirsions exhibited for unbalances of 0. 93 and 1. 67 gram -centimeters, 
tests at greater xinbalances were not attempted to avoid damaging rig components, such 
as the l^yrinth seals located at the bearing supports. 

Figure 9(b) is the filtered response of disk 3, located at the end opposite the turbine 
drive. Again, for the undamped run, there are three distinct amplitude buildups, which 
are at frequencies approximating those predicted by the tuidamped critical speeds pro- 
gram. The end disk tends to move out to amplitudes greater than indicated on this plot; 
however, the damper joturnal can move out only as far as the clearance in the bearii^ will 
permit, and this movement is approximately 127 micrometers. This may account for the 
jagged appearance of the trace at the second and third critical speeds. With the activa- 
tion of the dampers, the trace smooths out; again the theoretical prediction of the max- 
imum amplitude of 44 micrometers is low. 
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Low-speed measurements made at a number of places along the shaft from bearing to 
bearing indicated a bowing or warping of the shaft. The maximum amount of warp or bow 
at either side of the center disk was approximately 10 to 13 micrometers. This shaft 
runout can be seen in all the plots. The in-line, out-of -phase unbalance distribution was 
inadvertently placed in a plane only 40*^ away from the plane of the residual warp. This 
then could account for the shaft having a strong tendency to bend a greater amoimt in this 
plane than in any other. 
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Figure 10. - Comparison of rotor response for undamped and damped con- 
ditions. Type II unbalance distribution; unbalance, U, 0.62 gram- 
centimeter. 
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The rotor responses shown in figure 10 are for an in-line, alternating-phase distri- 
bution for two disks only; the center disk did not have an unbalance weight (see inset). 

This t 3 ^pe of unbalance distribution should have excited the second mode shape (see fig. 7). 
The undamped response, however, is similar to the previous in-line alternating phase 
response for three disks. Activating the oil-film dampers significantly reduced the res- 
onant amplitude. An amplitude of approximately 36 micrometers was recorded at 9700 
cycles per minute, whereas the imbalance response program predicted an amplitude of 
20. 3 micrometers at 10 000 rpm. 
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Figure 11. - Comparison of rotor response for undamped and damped con- 
ditions. Type III unbalance distribution; unbalance, U, 0.93 gram- 
centimeter. 



A total unbalance weight of 0. 154 gram is quite small compared with a rotor weight 
of 4. 54 kilograms. It is apparent from these tindamped rotor tests that this small 
amount of unbalance is not capable of inducing rotor bending in the classical mode shapes 
depicted in figure 7, particularly when both a residual bow and some residual imbalance 
(fig. 8) are present. 

The uiibalance distribution used to generate the ejcperimental responses shown in 
figures 11(a) and (b) for the center and end dif^, respectively, is an in-line, in-phase 
distribution using all three disks (see inset). The oU -squeeze -film dampers cause a re- 
duction in the peak-to-peak amplitude of over 80 percent. The predicted amplitude at the 
center disk of 30 micrometers at 9500 rpm for the damped rotor compares favorably 
with a measured amplitude of approximately the same magnitude at a rotor frequency of 
9600 cycles per minute. 

The response for an in-line, in-phase unbalance distribution applied to the test rotor 
correlated well with the predicted response at very small unbalances, as shown in fig- 
ure 11. It was decided to examine the response for much greater magnitudes of unbal- 
ance. In figure 12 the actual and theoretical rotor responses are compared. It is ap- 
parent in fig^ure 12(a), which is the peak-to-peak amplitude of the center disk, that 
squeeze-film dampers located at the bearings of a long, slender rotor are going to have 
little effect on attenuating midspan excursions. The agreement between the actual re- 
sponse and the predicted is quite good for both the amplitude and the frequency at which 
resonance occurs. 

The theoretical response program will predict, with considerably accuracy, the am- 
plitudes, phase angles, and forces transmitted if the input data, such as support flex- 
ibility and damping, are accurate. In this case the matchup is close. Stiffnesses of the 
centering springs were carefully measured, and this value was entered in the program 
as input. Damping coefficients were calculated from theory and for this case agreed 
quite well with measured values, at least at the resonant frequency. Damping values of 
6300 newton-second per meter shown on the plots were calculated for a concentrically 
positioned damper journal. 

It should be noted here that the recorded amplitudes of motion and transmitted 
forces were in the vertical or Y direction. Response was also measured in the X 
direction and appeared similar to Y measurements. The value of radial clearance used 
in calculating the damping for the theoretical unbalance response was 63. 5 micrometers, 
which was also the measured assembled clearance. 

In figures 12(b) and (c), considerable disagreement is apparent between the mea- 
sured amplitudes of motion and the theoretical for both the end disk and the damper jour- 
nal. The theoretical response is much greater than the measured response at frequen- 
cies above 10 000 cycles per minute. It is not exactly clear why such diagreement exists 
at the higher imbalances of 10. 6 and 15. 1 gram-centimeters. The rotor end and damper 
journal, opposite the turbine, deflect much less than predicted. The reason for this may 
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PEAK-TO-PEAK AMPLITUDE, (im 



ROTOR FREQUENCY, cpm 

(a) Center disk. (b) Disk 3. 


Figure 12. - Comparison of experimental and theoreUcal damped rotor response. Type III unbaiance distribution; unbalance. U, 10.6gram- 
centimeters; theoretical elastic force. K, 23.63x10^ newtons per meter; damping coefficient, 6300 newton-seconds per meter; radial 
clearance. Co. 63.5 micrometers. ^ 




be the noncylindrical shape of the damper bearing (see fig. 5). The walls of the bearings 
were thinner where distance probes were mounted. It is possible that this member was 
flexing asymmetrically, particularly with the higher rotating unbalance, and therefore 
the absolxite damper journal motions measured were less than if the damper bearing were 
perfectly rigid. The transmitted force also would be less in the actual system , since it 
is a function of the sum of elastic and damping forces, both of which, in turn, are func- 
tions of the displacement. 

The photographs of figure 13(a) show the damper journal orbital motion. The pic- 
tures on the left are for an undamped system, with the damper oil supply off. The pic- 
tures on the right are for identical conditions of unbalance, magnitude, and frequency ex- 
cept that oil was supplied to the dampers. It shoxold be pointed out that these pictvnes are 
of the unfiltered signals, and, therefore, Lissajous patterns of orbital motion are gen- 
erated which can be nonsynchronous as well as synchronous. The large orbit of the vin- 
damped journal motion at 9020 cycles per minute had a synchronous critical frequency ' 
with a peak-to-peak amplitude of 96.4 micrometers. The amplitude, frequency trace of 
figure 13(b) is a filtered signal made at the same time. The picture of the orbit was 
snapped while the rotor was approaching the second critical speed of 9600 rpm. The un- 
damped trace of figure 13(b) shows only resonances which are S 3 mchronous. Some of the 
multiple harmonics that occur in a lightly damped ball-bearing -supported rotor can have 
large amplitudes that can be quite destructive. These motions would not be seen at all if 
only the filtered trace were being observed. The Lissajous pattern at 25 710 rpm in fig- 
mre 13(a) has one inward loop; this signifies a miiltiple harmonic having a frequency twice 
that of the rotor speed. This orbit is a result of the gravity vector of a horizontally sup- 
ported rotor superimposed on the unbalance whirl vector (ref, 12). 

The noteworthy point to be made concerning these pictures is that oil-squeeze-film 
damping is very effective in attenuating the motions of both synchronous and nonsynchro- 
nous whirl. The amplitudes of motion for the damped tests were all less than those for 
the undamped tests. This fact is apparent in the filtered amplitude frequency traces 
shown in figure 13(b), The damped trace is superimposed on the undamped for a direct 
comparison of the synchronous motions. 


Oil -Squeeze-Film Damping 

Presented in figures 14 to 16 are oil-squeeze-film damping coefficients, the force 
transmitted through the oil film, and the eccentricity ratio, all plotted as functions of 
rotor frequency. Eccentricity ratio, as defined here, is the average of the peak joinnal 
displacements in the X-Y plane divided by the radial clearance in the damper. Table I 
contains a s\unmary of average damping coefficients at different magnitudes of unbalance 
and types of unbalance distributions. 




14 820 rpm 



18 600 rpm 


UNDAMPED RESPONSE DAMPED RESPONSE 

(a) Orbital motion. 

Figure 13. - Orbital motion of damper journal (bearing 2) with and without oil and corresponding 
amplitude, frequency response. Type I unbalance distribution; unbalance, U, 0.93 gram-centimeter; 
radial clearance. Of, 63.5 micrometer; oil flow, 5.15X10'^ cubic meter per minute; supply pressure, 
0.19 megapascal. 
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25 710 rpm 



29 500 rpm 


UNDAMPED RESPONSE 


(a) Concluded. 
Figure 13. - Continued. 


DAMPED RESPONSE 















The experimental values obtained in this investigation are averages calculated from 
amplitudes, forces, and phase angles measured in the X-Y plane. The relation of the 
force transducers and displacement sensors is shown in the schematic of figure 6. Be- 
cause displacement sensors and force transducers were not located in the same trans- 
verse planes, the phase angles measured had to be adjusted by 45°, This probably re- 
sulted in some error in calculating the dampii^ coefficients, but because the orbits were 
nearly circTilar (see fig. 13), the error introduced was not large. The straight lines 
plotted in figures 14 to 16 are the values of theoretical damping for a concentrically po- 
sitioned journal (eccentricity ratio 77 = 0) and a 180° film. 

The experimental values of transmitted force and eccentricity ratio plotted in these 
figures are also averages of values measured in the X-Y plane. These values are com- 
pared with those obtained in the unbalance response program. 

Values of damping coefficient were obtained by measurii^ the force transmitted 
through the oil film to the damper bearing F,pp^, the displacement of the damper journal 
e, the angle by which the displacement lags the exciting force (acts 180° from trans- 
mitted force) (p, the rotor frequency co, and the amplitude of the transmitted force Fq. 

The derivation of the equations used in the data reduction is given in the appendix; 
the equations are as follows: 

Damping coefficient: 


Stiffness coefficient: 


x,y “ 


e w 
x,y 


(Fp) cos (P 
K = V °/x,y 
x,y e 

x,y 


Stiffness coefficients are not presented in this report. The values calculated from 
this expression would depend on the cosine of the measimed phase angles. The accuracy 
of the angle near 90° is approximately ±1°. This error is not serious when the sine of 
the an^e is taken, since the sine does not vary significantly for 2° when near 90°. How- 
ever, the cosine does vary significantly even for 2°, and consequently there is a signifi- 
cant error introduced in calculated values of the oil -film stiffness coefficient. 
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Theoretical dampii^ coefficients plotted in figures 14 to 16 were calculated from the 
following eqTiation: 


B = 




o 2 

C|(l - r,^) 


where p is oil viscosity, R is journal radius, L is journal land length, Cj^ is radial 
clearance (diametral clearance/2), and rj is the conventional journal eccentricity ratio, 
defined as a static deflection Oq diAdded by the radial clearance Cj^. For the theoret- 
ical damping coefficients calculated in figures 14 to 16, rj is zero. Dynamic eccentricity 



ROTOR FREQUENCY, cpm 


(a) Unbalance, U. 0. 93 gram-centimeter. 

Figure 14. - Comparison of theoretical and experimental 
damping coefficients, transmitted forces, and eccen- 
tricity ratios as fu notions of rotor frequency. Type I 
unbalance distribution. 


25 



ratio 6 plotted in these figures is defined as the radius of the orbit e divided by the ra- 
dial clearance Cp. This equation is derived from the Reynolds lubrication equation for 
the case of a 180° film (cavitated bearing) by using short -bearing theory (ref. 10). 

Results plotted in figures 14(a) and (b) are for an in-line, alternating-phase unbal- 
ance distribution. The total unbalances were 0. 93 and 1. 67 gram -centimeter, respec- 
tively, Average values for the experimental dampii^ coefficient were 6740 and 6531 
newton -second per meter, as compared with values of 6492 and 6213 newton-second per 
meter for the theoretical, 180° film. In all the tests the damper oil supply pressure was 
set at 0. 19 megapascal, and the flow rate was 5. 15X10"^ cubic meter per min\ite. The 
fact that the experimental damping coefficients over the frequency range were slightly 
higher than the theoretical indicated that at certain conditions of speed and load the actual 
oil-film extent in the damper maybe slightly greater than 180°, 

The results plotted in figure 15 are for an in-line, alternating -phase unbalance dis- 
tribution for the two end disks only, designated unbalance type II. The experimental 
values of damping coefficients have somewhat more scatter than the previous data; how- 
ever, the values do not deviate too much from those for the theoretical 180° cavitated 
film predicted by short -bearing theory. Average values are compared in table I, Ex- 
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(b) Unbalance, U. 1.67 grani'centlmeters. 
Figure 14. - Concluded. 


26 



EXPERIMENTAL RESULTS 



8000 



o 

D 

A 

DAMPING COEFFICIENT 
TRANSMITTED FORCE 
ECCENTRICITY RATIO 
Tupnpv 

o z ■=& 

ooo 

iscPrlm 


1 ntUKT 

o O vt 

§ tt Z 

^ LU — 

<=» O - 

O CO 

6000 

4000 

I°1 

“o o 5“ 
1 1 1 ...J 

o 

LJ 

— o 'J 

1 1 1 1 1 1 



(a) Unbalance, U, 0.62 gram-centimeter. 

Figure 15. - Comparison of theoretical and experimental damping 
coefficients, transmitted forces, and eccentricity ratios as func- 
tions of rotor frequency. Type II unbalance distribution. 


perimental values of dampii^ coefficient deviate from theory most at rotor frequencies 
near the critical, which is approximately 9600 cycles per minute. Phase angles ap- 
proaching the critical frequency change very rapidly, and, therefore, the angles are sub- 
ject to more error than those angles measured at frequencies greater than the critical 
frequency. Rotor orbits at tiiese frequencies are also elliptical. 

Agreement between theory and experiment for the transmitted force and the eccen- 
tricity ratio is quite good throughout most of the frequency range. Experimental values 
tend to increase at a faster rate at higher rotor speeds. Ehridence indicates that a fourth 
bending critical is being approached. 

In figure 16 damping coefficient, transmitted force, and eccentricity ratio are 
plotted as functions of rotor frequency. The unbalance distribution used in the three 
disks was an in-line, in-phase distribution, designated type m. Magnitudes of total un- 
balance varied from 0. 93 to 15. 1 gram -centimeter. The experimental values vary from 
6. 8 percent below to 9. 2 percent above the predicted values. The trend, however, indi- 
cates experimentally determined damping coefficients for these series of tests are in 
excellent agreement with values calculated from short -bearing theory for a 180°, cav- 
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(b) Unbalance, U. 1.11 gram-centimeters. 
Figure 15, - Continued. 


itated film. Experiments also indicate that damping coefficients are not a strong function 
of dynamic eccentricity ratio but, in general, remain fairly constant. Similar results 
are presented in reference 13; that is, squeeze-film damping coefficients, although 
greater than predicted by 180° film theory, remained fairly constant at dynamic eccen- 
tricity ratios varying from 0. 03 to 0. 4 for an end fed damper with a length -diameter ratio 
of 0. 34. 

The data plotted in figures 16(d), (e), and (f) are for total unbalances of 4. 97, 10. 6, 
and 15. 1 gram -centimeters, respectively. The resulting higher dsmamic eccentricity 
ratios were expected to result in higher damping coefficients. This did not occur; in 
fact, an opposite trend was observed, in that damping coefficients actually decreased. 
Also, a hydrodynamic instability was induced. This type of operation is referred to in 
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ROTOR FREQUENCY, cpm 
(c) Unbalance. U. 1. 5 gram-centimeters. 

Figure 15. - Continued. 

the literature (refs. 14 and 15) as bistable operation. Briefly, it is operation of the 
damper journal in a stable condition at either of two eccentricity ratios. The region of 
frequency where this instability occurred is shown in figures 16(e) and (f). The rapid in- 
crease in journal amplitude is also shown at approximately 25 000 rpm in figure 12(c). 

A region is indicated since the bistable phenomenon noted did not always occur at exactly 
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ROTOR FREQUENCY, cpm 
Id) Unbalance, U, 3. 3 gram-centimeters'. 

Figure 15, - Concluded. 

the same frequency in subsequent trials. It is not the purpose of this investigation to 
examine non-steady-state behavior other than to note its dramatic effect in the reduction 
of the damping coefficient. The oil flow to the damper increased at frequencies approach- 
ing the bistable region. It is postulated that the damper at this juncture became starved, 
and the oil film ruptured. This would account for the decrease in damping coefficients 
above the instability region. 
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(a) Unbalance, U. 0.93 gram-centimeter. 

Figure 16. - Comparison of theoretical and experimental damp- 
ing coefficients, transmitted forces, and eccentricity ratios 
as functions of rotor frequency. Type III unbalance dis- 
tribution. 
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SUMMARY OF RESULTS 


Ejqierimentally determined values of squeeze-film dampii^ coefficients were in good 
agreement with theoretical values calculated from equations derived by the short -bearing 
approximation for a 180° film and a concentrically positioned journal. Values remained 
fairly constant at rotor frequencies up to 31 000 cycles per minute. 

Central -feed-groove oil-squeeze-film dampers, located at the bearings, were most 
effective in reducing amplitudes of motion at the rotor ends. They were, however, less 
effective in attenuating amplitudes of motion at the rotor center. A reduction in ampli- 
tude of 88 percent was observed for an imbalance moment of 0. 93 gram -centimeter dis- 
tributed in line and in phase over three disks. 

When oil-squeeze-film dampers were activated, for an in-line, alternating-phase im- 
balance distribution, three distinct resonant amplitudes were reduced to only one. The 
mode shape for this one corresponded to the third bending critical speed. 

The oil -squeeze -film dampers were effective in reducing the amplitudes of motion 
for nonsynchronous whirl as well as synchronous. 

An unbalance-response computer program was foimd to be quite accurate in predict- 
ing the response of the three-disk ejqperimental rotor used in this investigation. The dis- 
crepancies that did occur apparently were due to the presence of a residual bend in the 
rotor havii^ a shape similar to that which occurs at the first critical speed. The place- 
ment of weights in a line and near the plane of residual bend induced greater bending than 
would probably have occurred with a straighter shaft. Conversely, locating the center 
weight 180° from the weights in the two end disks resulted in the rotor not bending as 
much as predicted by theory. 

Unbalances of 10. 6 and 15. 1 gram -centimeters arranged in line and in phase in all 
three disks produced an instability in the rotor referred to as bistable operation. Damp- 
irg coefficients were reduced significantly at frequencies greater than the onset fre- 
quency. Damper oil flow rates increased, probably as a result of the higher damper 
journal eccentricities and increased hydrodynamic pumping. 

With the undamped rotor, resonant amplitudes occurred at frequencies slightly 
greater than those predicted by a critical-speed, mode-shape program. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, Jtily 28, 1977, 

505-04. 


APPENDIX - DERIVATION OF DATA REDUCTION EQUATIONS 

The force transmitted through the oil film to the damper bearing, and measured by 
the force transducers, is equal to the sum of the elastic force, the dissipative or damp- 
ing force, and the oil-film inertial force: 

transmitted force = elastic force + dampir^ force + inertial force 

Because the volume of oil in the squeeze film is so small, the inertial force can be 
neglected, and 





d(x,y) 

dt 


(Al) 


It is assumed that the journal processes about the bearing center in an elliptical orbit and 
the motion is harmonic: 


X = e^ cos a»t 


y = Cy sin cot 


— = w sin cot 
dt ^ 

^ = e„ CO cos cot 

dt y 




(A2) 


(A3) 


Let (Pq) be the amplitudes of the transmitted force in the x and y directions, re- 
V vx,y 

spectively; then 



Fy = sin(cot + <?y) J 

Substituting ecjuations (A2), (A3), and (A4) in (Al), we obtain 
Fj^ = ^ ~ ®x®x'^ ^ ~ (^o) 


(A4) 


(A5) 
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and 


Fy = Ky6y sin cot + ByCyCo cos cot = ^FqJ sin(cot + <Py) 


(A6) 


where (p and <p are the angles by which the displacement lags the exciting force (see 
X y 

sketch (a)). 


Y 



Let y = 0, X = e , and cot = 0; solvii^ for the stiffness and damping coefficients gives 


= (^ol fx' 


J^x = 




COS (p. 


X 


(A7) 


Ky6y(0) + ByeyCO(l) ^ Sm(0 + <j9y) 


y 

sin w 


y eyco 


(A8) 


Let X = 0, y = ®yj and cot = tt/ 2; again solving for the stiffness and damping coefficients 
gives 


Vx<»> - Vx“<l' = (I + 1’x) 


®x = 




sin (p 


(A9) 


e^co 
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Illlllllllllllllllllll I 


K(ey)(l) +BeyW(0) = (Fo)^ 


cos «Py 



(AlO) 


The transmitted forces in the X-Y plane are measured along with the journal dis- 
placements in orthogonal planes. The rotor frequency o) is also measured. When these 
values are known, the dynamic film coefficients in the X-Y planes can be solved for. 
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